o 

O 
in 



o 
o 



On 
ON 



X 



Dispersionless Fermionic KdV 



J. Barcelos-Neto 
Instituto de Fisica, 
Q\ [ Universidade Federal do Rio de Janeiro, 

Rio de Janeiro, RJ 21945-970 
Brazil 



and 



Alin Constandache and Ashok Das 
Department of Physics and Astronomy, 
O University of Rochester, 

Rochester, New York, NY 14627-0171 
USA. 



Abstract 



o 

We analyze the dispersionless limits of the Kupershmidt equation, the SUSY KdV-B equation 
and the SUSY KdV equation. We present the Lax description for each of these models and bring 
out various properties associated with them as well as discuss open questions that need to be 
addressed in connection with these models. 



1 Introduction: 



In the last two decades, integrable models [1-3] have drawn a lot of attention from various points 
of view. More recently, however, research in a particular class of integrable models, known as 
the dispersionless limits of integrable models, has become quite active. These involve equations 
of hydrodynamic type [4-8] and include systems like the Riemann equation [8-10], the polytropic 
gas dynamics [8,11], the chaplygin gas, the Born-Infeld model [12-13] etc. These are models which 
can be obtained from a "classical" limit [6,14] of integrable models where the dispersive terms are 
absent. Thus, the dispersionlesss limit of the KdV equation [15], for example, corresponds to the 
Riemann equation and so on. The interesting thing about these models (and the "classical" limit) 
is that, if we know the Lax operator description of a given model, then, the dispersionless limit 
also has a Lax description, in terms of a Lax function on the phase space, where the commutator is 
replaced by a Poisson bracket. The Lax function can, in fact, be determined from the Lax operator 
in a systematic manner [7]. (Of course, there are also dispersionless models whose integrable, 
dispersive counterparts are not known.) 

While a lot is known about the dispersionless limits of bosonic theories, nothing is known 
about the fermionic theories (with or without supersymmetry). The main difficulty lies in the fact 
that, in dealing with fermionic models, we have to deal with a phase space with both bosonic and 
fermionic coordinates. The classical fermionic variables are nilpotent and, consequently, a power 
series representation that is so crucial in a Lax description seems to fail. Thus, for example, the 
supersymmetric KdV equation [16] is described in terms of a Lax operator which involves powers of 
the supercovariant derivative D. This is a fermionic operator whose square equals d (= ^). On the 
other hand, a phase space description of the model in the dispersionless limit would naively seem 
to require using a fermionic phase space variable which would give the same action as D through a 
Poisson bracket relation. Such an object, however, would be nilpotent and there cannot be a power 
series in this variable. Such difficulties have hindered the understanding of such fermionic theories 
so far. 

In this paper, we make the first attempt towards understanding such models. We analyze the 
dispersionless limits of fermionic KdV equations which include the Kupershmidt equation [17], the 
SUSY KdV-B equation [16,18] as well as the SUSY KdV equation [16]. All these equations can 
be thought of as fermionic extensions (with or without supersymmetry) of the Riemann equation 
which are integrable (Of course, we only concern ourselves with N = 1 supersymmetry). We obtain 
a "classical" Lax description for each of these models. However, the Lax functions and the Lax 
equations are obtained by brute force (with a lot of hard work) and we do not yet understand a 
systematic way in which the phase space Lax functions (for the fermionic models) can be constructed 
from the corresponding Lax operators. This remains an open question. From the Lax description, 
we construct all the local conserved quantities for the models in the standard manner. However, 
as we emphasize in the paper, fermionic models, in particular, the supersymmetric ones contain 
nonlocal conserved charges [19-20] as well, and it is not at all clear how a phase space Lax description 
can generate such quantities from a "Trace" of the Lax function. This, too, is an interesting open 
question. While the Lax description of the bosonic models in the dispersionless limit does give 



the first Hamiltonian structure from a generalized Gelfand-Dikii bracket [10], as we discuss, except 
for the SUSY KdV-B equation, we are unable to obtain the first Hamiltonian structure from a 
generalized Gelfand-Dikii bracket for the other two models. This is in spite of a Dirac analysis 
which we have carried out for the constrained form of the Lax functions [21]. In the case of SUSY 
KdV, this may be understood as signifying that the first Hamiltonian structure of SUSY KdV 
vanishes in this "classical" limit. However, for the dispersionless Kupershmidt equation, there does 
exist a first Hamiltonian structure and it is unclear how to obtain this from the Lax description 
itself. Our paper is organized as follows. In section 2, we recall, very briefly, the essential features of 
the dispersionless KdV equation. In section 3, we present the Lax description for the dispersionless 
Kupershmidt equation and discuss all its properties in detail. In section 4, we discuss the analogous 
results for the dispersionless SUSY KdV-B equation. Here, there are nonlocal conserved charges as 
well, and we present the algebra of the charges which takes a particularly simple form. In section 
5, we discuss results for the dispersionless SUSY KdV equation. Here, too, there are nonlocal 
charges and we present the algebra of the charges in a closed form. In section 6, we present a short 
conclusion emphasizing the open questions within the context of such models. 



2 Dispersionless Limit of KdV Equation: 



In this section, we will briefly review the known features of the dispersionless KdV equation [8,10]. 
As is well known, the KdV equation [15] 



du du d 3 u . . 

dt dx dx 3 



can be described by the Lax equation 



where the Lax operator 



§ = 4[(L 3 / 2 ) + ,L] (2) 

L = d 2 + u{x) (3) 

with d representing ■£-, u the KdV variable and ()+ denoting the differential part of a pseudo- 
differential operator. The third derivative term on the right hand side of eq. (1) represents a 
dispersive term and the dispersionless limit of this equation (namely, the equation where the dis- 
persive term is absent) is obtained as follows. Let 

d d d d 

dt dt' dx dx 

without rescaling the dynamical variable. Then, in the limit, e — > 0, the KdV equation reduces to 



which is the Riemann equation and has no dispersive term. This is known as the dispersionless 
limit of the KdV equation and it has a Lax description as well. Consider the Lax function on the 
classical phase space 

L=p 2 + u(x) (6) 

where p is the momentum variable of the phase space. Then, with the canonical Poisson bracket 
relations, it is easy to check that 

§ = -4{(L^) + , L } (7) 

with ()_|_ now representing terms with non-negative powers of p, gives the Riemann equation or 
the dispersionless limit of the KdV equation. Thus, one can think of the dispersionless limit 
as sort of a "classical" limit where the Lax operator goes into a Lax function {d — > p) and an 
operator Lax equation involving a commutator goes into a Lax equation involving Poisson brackets. 
Furthermore, the conserved quantities and the Hamiltonian structures (at least, the first structure) 
can be obtained from this Lax function [10]. (The difference in sign between the operator equation 
(2) and the phase space equation (7) results from the fact that [d,f] = 4- while {p,f} = — 4z- 
This would be reflected in all the cases that we discuss.) 



3 Dispersionless Limit of Kupershmidt Equation: 



The Kupershmidt equation [17] is a nontrivial fermionic extension of the KdV equation which is 
integrable. The dynamical equations, in this case, involve a bosonic variable u as well as a fermionic 
variable ip and are given by the coupled equations 

iH = Uxxx + 6uu x - 12ipip xx 

ipt = ^tpxxx + 3u x il> + 6uip x (8) 

where the subscripts denote differentiation with respect to the particular variable. This is not a 
supersymmertric system, but gives a nontrivial, coupled boson-fermion system (reducing to the 
KdV equation when ip = 0) which has a bi-Hamiltonian structure and is integrable. 

The Kupershmidt equation can be described by a Lax equation. Consider a Lax operator of 
the form 

L = d 2 + u + i}d- 1 i> (9) 

Unlike the KdV equation, this Lax operator is truly a pseudo-differential operator involving the 
fermionic variable ip. Once again, it can be easily checked that the Lax operator equation 

§ = 4[(L 3 / 2 ) + ,L] (10) 

gives the Kupershmidt equations. 



In trying to derive the dispersionless limit of the Kupershmidt equation, one runs into various 
problems. First, we note that, under the scaling d — ► ed, 

d- x ^e- l d (11) 

so that in the limit e — > 0, the term in the Lax operator containing the fermionic variables would 
appear to diverge. Second, if we naively let d — > p, then, of course, the fermionic term in the 
Lax operator would vanish. An alternate approach is to recognize that d~ l is really an operator 
which can be taken to the right with the help of the Leibnitz rule, giving an infinite series of 
terms in which one can let d — > p and the terms involving the fermionic variables would no longer 
vanish. However, a short calculation shows that such a procedure leads to an inconsistent Lax 
equation. Therefore, finding a Lax function and a Lax description for the dispersionless limit of 
the Kupershmidt equation genuinely poses a challenge. 

The solution to this problem comes as follows. Consider the Lax function 

L = p 2 + u-p~ 2 ipip x (12) 

Namely, the term containing the fermionic variables corresponds only to the first nontrivial term 
in taking the operator <9 _1 to the right and setting d = p. It is now easy to check that the Lax 
equation (once again, note the difference in sign) 

§ = -4{(L 3 / 2 ) + ,L} (13) 

gives rise to the equations 

u t = 6uu x — 12ipip xx 

ip t = 3u x ip + 6uip x (14) 

which indeed represent the dispersionless limit of the Kupershmidt equation in (8). (It can be 
checked that the "classical" limit involves the scaling d — > ed and tp — > e" 1 ' 2 tp without which 
the fermion terms would not be present in the boson equation. This is indeed very different from 
the clasical limit in bosonic theories. We will see this in all the fermionic models we describe.) 
The passage from the Lax operator description of the Kupershmidt equation to the Lax function 
description of the corresponding dispersionless equation clearly is not as straightforward as the 
bosonic KdV equation and it is not at all clear how the Lax function of eq. (12) could have been 
deduced from that of eq. (9). 

The existence of a Lax description, as we know, gives, in a simple way, many of the interesting 
properties of the integrable system. Thus, for example, from the Lax function, we can define 



H n = TrL— 



2n+i . ...u-i / — — \ / . n+ i , i\„,n-1 




u n+i - 2n(n + l)u n - L ^ x (15) 



where "Tr" represents the integral of the coefficient of the p 1 term in the expression. We can, in 
fact, easily check, using the equations of motion, that the i? n 's are conserved. Furthermore, these 




can be identified with the limit e — ► of the conserved charges for the Kupershmidt equation under 
the scaling 

d -» ed, V -> £~ 1/2 4> (16) 

This is again a manifestation of the necessity for scaling the fermion variables to obtain the disper- 
sionless limit (without the scaling, there would be no fermion terms). 

The Lax description, for the bosonic KdV equation, of course, gives the Hamiltonian structures 
(at least the first one) naturally through a generalization of the Gelfand-Dikii bracket [10]. A similar 
analysis fails in this case, in spite of a careful treatment (Dirac analysis) of the constrained nature 
of the Lax function [21-22]. The derivation of the Hamiltonian structure for the Kupershmidt 
equation from the Lax description, therefore, remains an open question. However, the Hamiltonian 
structures for this system are not hard to obtain directly. In fact, it can be checked easily that 

V 1 = 

give the first two Hamiltonian structures of the system which can also be derived from the first two 
Hamiltonian structures of the Kupershmidt equation under appropriate scaling (see (16)). Once 
we know the Hamiltonian structures, it is straightforward to show that the conserved quantities 
are in involution. For example, with the first Hamiltonian structure in (17), we have 

{H n ,H m } 1 = -2m(n+l){m+l)(m-l)Jdx(u n+m - 2 ^ x ) x 

= (18) 

with the usual assumptions on asymptotic fall off of the dynamical variables. This proves that the 
model remains integrable in the dispersionless limit. 



4 Dispersionless Limit of SUSY KdV-B Equation: 



The super symmetric KdV-B equations correspond to a trivial supersymmetrization of the KdV 
equation [16,18] and are given by 

ut = Uxxx + 6uu x 

4>t = ?Pxxx + Quip x (19) 

This is a set of simple equations where the boson equation does not depend on the fermionic variable. 
However, the pair of equations in eq. (19) are invariant under the supersymmetry transformations 

St/} = Xu 

5u = -\il> x (20) 



with A a constant Grassmann parameter of the transformation. It is this particular supersym- 
metrization of the KdV equation which shows up in the supersymmetric one matrix model. 

The SUSY KdV-B has a Lax description and, being a supersymmetric system, the proper 
description for it is in superspace. Let us define a fermionic superfield 

$(x,0) = ij){x) + 6u{x) (21) 

where 9 represents the Grassmann coordinate of the superspace. Let us now define a Lax operator 

L = £> 4 + (D$) (22) 

where the supercovariant derivative is defined as 

D -w +e h D '' = B (23) 

It is now straightforward to check that the Lax equation 

f t =A[(L^) + ,L] (24) 

where () + refers to non-negative powers of D, leads to 

$t = $xxx + HD$)<S> X (25) 

which contains both bosonic and fermionic components of the SUSY KdV-B equation. 

This Lax operator is in many ways reminiscent of the Lax operator for the KdV equation and 
finding the dispersionless limit of this system, therefore, does not pose too much difficulty. However, 
there are some interesting features that arise in the case of the dispersionless equation which we 
will describe. First, let us note that although the Lax operator for the SUSY KdV-B system is 
written in terms of the supercovariant derivative D, we could have written it equally well in terms 
of d because of the relation between the two. Keeping this in mind, let us look at the Lax function 

L = p 2 + (D$) (26) 

It is then, easy to check that the Lax equation 

f t =-A{{L^) + ,L] (27) 

gives the dispersionless limit of the SUSY KdV-B equation, namely, 

*t = 6(£>$)$ x (28) 



which, in components, takes the form 



ut = 6uu x 

4>t = 6u^ (29) 



These equations are indeed the dispersionless limits of the SUSY KdV-B equations in (19) and we 
can think of them as the trivial supersymmetrization of the Riemann equation. 



The Lax description, as we have seen, gives us the conserved quantities and it is easy to check 



that 



H„ 



n+l 



2 n + 2 U\ 2n+3 

sTr L 2 



(2n + 3)!! 
1 



(n + l)(n + 2) 
! f dxu n+1 ^ 



dz {D$) 



ra+2 



n + l 



(30) 



are conserved for n = 0, 1, 2, ■ ■ ■. Here "sTr" stands for the integration of the coefficient of p~ l over 
the entire superspace with dz = dxdO. We have also chosen a particular normalization for later 
purposes. There are several things to note about these conserved charges. First, since they are 
expressed as integrals over the superspace, they are automatically invariant under supersymmetry 
transformations. However, we see that these charges are fermionic in nature and this suggests that 
the Hamiltonian structure for the system should be odd (anti-bracket structure). This can, in fact, 
be readily verified, namely, let us define the dual 



Q = <?o + qip 1 



so that 



(31) 

(32) 



sTrLQ= f dz qi (D<£>) = - f dz (D qi )$ 
Then, a generalization of the Gelfand-Dikii definition of the first Hamiltonian structure [10] 

{sTr LQ, sTr LV] = sTr L {Q, V} (33) 

yields 



{$(z), $(/)} = -25{z-z') 



(34) 



which is, indeed, the correct Hamiltonian structure for the system (up to normalization) and is an 
odd structure. This is a special feature of the SUSY KdV-B system [23] (as well as the SUSY TB-B 
[24] system). 



One way to understand the odd Hamiltonian structure is as follows. Let us consider the La- 
grangian 



L 



dz 

dx 



1 



2 (n + l)(n + 2) 



(Z><&) 



n+2 



1 



Ulpt ~ ipu t ) 



1 



n + l 



n n+i V, 



(35) 



This Lagrangian describes the dispersionless limit of the SUSY KdV-B hierarchy as its Euler- 
Lagrange equations and is clearly fermionic. It is then, easy to work out from the special structure 




of this Lagrangian that there are constraints in such a theory which modify the canonical even 
Poisson brackets to Dirac brackets which are odd. 

The first Hamiltonian structure can be written as a 2 x 2 matrix in terms of the components 
as (up to a normalization) 

' ° -' ) (36) 

It is also easy to check from the explicit forms of the conserved quantities in eq. (30) that the 
recursion operator (function) for the system is given by 

so that 

/ SH n+2 \ / SH n+1 \ 

afo, )=R( «fcx (38) 

V Si> ) \ Sip J 

The recursion operator, then, leads to the hierarchy of Hamiltonian structures for the system given 
by 

Vn+l = R n V l= ^ n ~ Un \ (39) 

and all of them are indeed odd structures. 

In addition to the fermionic conserved charges, H n+ \, the dispersionless SUSY KdV-B equation 
also has bosonic conserved charges. It is easy to check that the charges 

Q n+l = jdz$(D$) n 

dx(u n+1 -nu n ~ l ^ x ) (40) 



are conserved. Furthermore, these charges are supersymmetric and can be obtained from the 
conserved bosonic charges of the SUSY KdV-B system (which are obtained by taking "super Trace" 
of multiples of quartic powers of the Lax operator). However, the charges, in the present case, cannot 
be obtained from the Lax function of the system. This is, in fact, a puzzling general feature of such 
systems. Namely, we know that supersymmetric integrable systems possess nonlocal conserved 
charges and they arise, say in the case of SUSY KdV, from the trace of quartic powers of the Lax 

2n + l 

operator such as sTr L 4 . In the case of operators, this is meaningful because the supercovariant 
derivative, D, is an operator square root of d. In a classical phase space, however, no such relation 
exists and it is not clear how to obtain the analog of the nonlocal charges in such a case. Thus, for 
example, in the dispersionless SUSY KdV-B system, it is easy to see that the quantities 



dxu n = I dzD-\D<5>) n (41) 

are conserved. However, these are nonlocal quantities and it is not clear how to obtain such 
conserved quantities from the Lax function. Let us note that g n 's are purely bosonic, and are not 



invariant under super symmetry transformations (even though they can be written as an integral 
over superspace, they involve D^ 1 which leads to violation of supersymmetry). 

Finally, let us note that we have a set of bosonic and fermionic conserved charges, in this theory, 
which are invariant under supersymmetry. Therefore, it is meaningful to investigate the algebra 
of these charges. With the first Hamiltonian structure, for example, we can easily calculate the 
algebra which has the form 

{H n ,H m } = 
{Qn,H m } = 
{Q n ,Qm} = (n-m)(n-l)(m-l)fl n+m _i (42) 

The first equation simply says that the conserved Hamiltonians are in involution so that the system 
is integrable. The second is an expression of the fact that the Q n 's are conserved under all the 
flows of the hierarchy. The nontrivial and really interesting one is the last one which is reminiscent 
of supersymmetry algebras. The important thing to remember in this algebra is the fact that Q n 's 
are bosonic, H n 's fermionic and the bracket is odd. 



5 Dispersionless Limit of SUSY KdV Equation: 



The super symmetric KdV-B equation represents a trivial supersymmetrization of the KdV equation 
and hence there was not much difficulty in taking its dispersionless limit. The N = 1 supersym- 
metric KdV equation, on the other hand, is a case where we expect some challenge in taking the 
dispersionless limit, just like the Kupershmidt equation. Let us recall that the equation [16] 

$ t = (D 6 <S>) + 3D 2 ($(D<5>)) (43) 

which in components has the form 

u t = u xxx + 6uu x - 3ipvp xx 

ipt = i>xxx + 3(uip) x (44) 

represents the N = 1 supersymmetric KdV equation which is known to be integrable. This equation 
is invariant under the supersymmetry transformations of eq. (20) and has a Lax description of the 
following form. Let us consider a Lax operator of the form 

L = D 4 + D$ (45) 

Then, it is easy to check that (in this case, there is a degeneracy and the Lax operator L = D 4 — &D 
works equally well) the Lax equation 

d ± = A[(L^) +1 L] (46) 
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gives us the SUSY KdV equation. It is here that we see the fermion nature of the problem coming 
into play. Unlike the case of SUSY KdV-B equation, here the Lax operator cannot be written 
completely in terms of the bosonic d. One has to understand how to take the "classical" limit 
of the operator D. Of course, classically, we would have a phase space which would involve both 
bosonic and fermionic coordinates. If p and p$ represent the bosonic and the fermionic momenta 
respectively, then, we can define a fermionic function 

U = -(p e + 0p) (47) 

which would generate covariant differentiation through [25] 

{U,A} = (DA) (48) 

for any arbitrary superfield A, and would also satisfy 

{U,U} = -2p (49) 

reminiscent of the operator relation (23). However, II is a fermionic function and, therefore, nilpo- 
tent, which is not quite the behavior of the operator D. 

Furthermore, if we look at the supersymmetry transformations of eq. (20), it is clear that the 
scaling d — > ed would lead to (in the limit e — > 0) 

5ip = An, Su = (50) 

which gives rise to a nilpotent algebra. This is also clear from an analysis of the supersymmetric 
algebra, namely, if it is only the generator of spatial translation which scales, then, the supersym- 
metry algebra would become nilpotent. To have supersymmetry as we know it, in the "classical" 
limit, we need, therefore, to scale the Grassmann coordinates and the fermionic variables as 

^ e -l/20 5 ^ _+ e -l/2^ ( 51 ) 

This is another indication that, unlike the bosonic variables, fermion variables need to scale for a 
consistent "classical" limit. 

From the structure of the Lax operator in eq. (45), we note that if we introduce, in analogy 
with the KdV equation, the Lax function 

L = p 2 - $11 (52) 

where II is defined in eq. (47), it is easy to verify that the Lax equation leads to inconsistencies. 
Thus, a priori, it is not clear how to proceed in finding a Lax description and the dispersionless 
limit. A little bit of analysis shows that the structure of the Lax function, in this case, is likely to 
be an infinite series of the form 

oo 

L = p 2 + J2( A n + nB n )p- n (53) 

n=0 

11 



where A n and B n are superfields to be determined recursively. An infinite series with undetermined 
coefficients, of course, would not be very useful in studying the properties of the system. 

Surprisingly, however, we have found that the finite order Lax function 

1 ~ 2 —4 

L = p 2 + - (£>$) + ^ ((D®) 2 - 2$3 X ) - V — $(£>$)$, (54) 

leads to the dispersionless limit of the SUSY KdV equation 

$ t = 3-D 2 ($(£>$)) (55) 

through the Lax equation 

§ = -4{(L 3 / 2 ) + ,L} (56) 

This is quite a nontrivial Lax function and it is interesting to note that it does not involve the 
fermionic variable IL It is described completely in powers of p, and it is not clear how one could 
have deduced this particular form of the Lax function from the Lax operator of SUSY KdV in eq. 
(45). 

On the other hand, since we have a Lax description of the dispersionless SUSY KdV, we can 
determine the conserved quantities of this theory. It is easy to check that the quantities 

2 +i 1 / 2n+l \ r 

H n = TrL^ = -[ 2 dz<&{D<S>) n (57) 



are conserved under the evolution of the system. It is interesting to note that, although the con- 
served charges are obtained as traces, they can actually be written in a completely supersymmetric 
manner and they are bosonic unlike the case in SUSY KdV-B. Unfortunately, the Gelfand-Dikii 
brackets do not give a meaningful first Hamiltonian structure for this system. While this may seem 
reminiscent of the Kupershmidt equation, in this case, it is easy to check that the first Hamiltonian 
structure of SUSY KdV vanishes under the scaling (with e = 0), namely, in the "classical" limit. 
On the other hand, it is not hard to obtain a Hamiltonian structure of the theory directly (as far 
as we know, it is only the first Hamiltonian structure which is known so far to be derivable from 
the Gelfand-Dikii brackets in the dispersionless limit of bosonic theories) and has the structure 

T> = -- (3$D 2 + (D&)D + 2(L> 2 $)) 5{z - z) (58) 

This, indeed, coincides with the "classical" limit of the second Hamiltonian structure of SUSY 
KdV. 

Once we have the Hamiltonian structure, it is easy to check that 

{ H n ,H m} = J S -§V S J^ = (59) 
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Namely, all the conserved charges are in involution as we would expect for an integrable system. 
In addition to these local conserved quantities, the dispersionless SUSY KdV also has two sets of 
nonlocal conserved quantities, namely, one can check directly that 

F n = J dz(D- l <b) n 

G n = I dz^D- 1 ®)* 1 (60) 

are conserved under the evolution. It is worth noting that while F n 's are fermionic, G n 's are all 
bosonic and that 

G = 4H (61) 

Once again, it is not clear how one can derive such nonlocal charges from a Lax function description. 
The algebra of these charges is tedious, but can be calculated with the Hamiltonian structure of 
the system in eq. (58) and has the form 

{H n ,H m } = = {F n ,H m } = {G n ,H m } 
\F n ,G m } = = {G n , G m \ 
{F n ,F m } = nmG n+m -2 (62) 

Vanishing of the first three brackets is simply a statement that the Hamiltonians, H n are in in- 
volution and that F n and G n are conserved. However, from the other relations, we see that F n 
and G n together define a supersymmetry algebra. In particular, it is not difficult to check that 
F\ corresponds to the generator of supersymmetry of the system (up to normalization) with Go 
(= 4:Hq) representing the Hamiltonian. 



6 Conclusion 



Although a lot is known about the structure of dispersionless limits of bosonic integrable systems, 
the dispersionless limits of fermionic integrable models have proven difficult for a variety of reasons. 
In this paper, we have made the first attempt towards understanding such models. We have studied 
the dispersionless limits of the Kupershmidt equation, the SUSY KdV-B equation as well as the 
SUSY KdV equation all of which can be thought of as fermionic extensions (with or without 
supersymmetry) of the Riemann equation. We have obtained the Lax function description for each 
of these systems and constructed the Hamiltonians for such systems from them. However, there 
is as yet no systematic understanding of how such Lax functions can be deduced from the Lax 
operator description of the original theory. We have presented the Hamiltonian structures, the 
nonlocal conserved charges for such systems, as well as the algebra of conserved charges in closed 
form. However, it is not clear how a Lax function description can give rise to nonlocal conserved 
quantities nor is a generalization of the Gelfand-Dikii brackets for the Kupershmidt and SUSY KdV 
equations is known. While we have discussed many interesting features of such fermionic models, 
we have also brought out several interesting open questions that deserve further study. 
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